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ABSTRACT 

We report high resolution (R^37,000) integral field spectroscopy of the central region (r<14arcsec) 
of the Red Rectangle nebula surrounding HD 44179. The observations focus on the A5800A emission 
feature, the bluest of the yellow/red emission bands in the Red Rectangle. We propose that the 
emission feature, widely believed to be a molecular emission band, is not a molecular rotation contour, 
but a vibrational contour caused by overlapping sequence bands from a molecule with an extended 
chromophore. We model the feature as arising in a Polycyclic Aromatic Hydrocarbon (PAH) with 
45-100 carbon atoms. 

Subject headings: planetary nebulae: individual (Red Rectangle), ISM: molecules, ISM: lines and bands 



1. INTRODUCTION 

The Diffuse Interstellar Bands (DIBs), first identified 
as interstellar absorption features by Merri ll and cowork- 
ers dMcrrill 1934, 1936; Merrill & Wilsoi] [1931 . repre- 
sent one of the longest outstanding mysteries in astro- 
nomical spectroscopy. 

The DIBs are widely believed to represent the origin 
bands of molecular transitions. Ultra-high resolution 
spectroscopy of a selectio n of the stronger DIBs con- 
firms their molecular origin l|Sarre et a L'iggSal lKerr etal\ 
(1998), while the lack of strong correlations between 
known bands precludes the identification of any tran- 
sitions to higher vibronic states. 

Identification of the carriers of the DIBs will ultimately 
require accurate measurement of the gas phase spec- 
trum of a range of extraterrestrial molecules to allow 
direct comparison with astro nomical observations (fo r 
a review of current work see ISchmidt fc SharnI ll200l : 
IFulara fc K rclowski (2003); IHerbid II1995D ). However, 
the parameter space occupied by the range of candidate 
carriers is prohibitively large for full investigation with 
the often complex techniques required to record the gas 
phase spectra. 

In order to restrict the parameter space one can turn 
to studies of alternate environments to the diffuse clouds 
within which most DIBs are observed. These alternate 
environments may represent the Rosetta stone for DIB 
spectroscopy-'^, offering the possibility of presenting ad- 
ditional features of the molecular species, providing the 
much needed additional constraints to the properties of 
the carrier molecules. 

In this paper we present a high spatial, and spectral, 
resolution study of one of the emission features seen in 
the Red Rectangle nebula. At a di stance of only 330pc 
(but see iMen'shchikov etoD ()2002j) for details of a re- 
Electronic address: rgsOaao. gov.au ^^^^^^^^ 
Electronic address: I'.bichmidtOchem. u syd.edu.au | 

^ We thank A. Witt for suggesting this elegant metaphor, private 
communication 



vised distance estimate of 710pc^) the nebula represents 
an excellent candidate for the study of the injection of 
material into the interstellar medium by a dying star. 

2. THE RED RECTANGLE 

The bi-conical Red Rectangle (RR) nebula 
l|Cohen etal\ 119751 l2004^ surrounding the evolved, 
mass-losing star HD 44179 is one of the strongest 
sources in the sky when observed in the mid-IR 
(IRAS 06176-1036), where a set of emission fea- 
tures commonly attribut ed to polycyclic aromatic 
hydrocarbons (PAHs, AUamandola etai\ 1)19851) 1 
dominate the spectrum. PAH molecules provide an 
elegant solution to a numbe r of astronomical issues 
( Uy exctinction curv e - dP onn. Hodge fc MentallJ 
1968); DIBs carriers - ll van der Zwet fc Allamand olal 
19851 t^egcr fc Dhcndcc mirt I I1985D:" "~UI B 
carriers IC^balle. Tielens. AUamandola et dA 119891) 1. 
With the coming of age of the new generation of 
mid IR instrumentation, PAH emission may prove to 
be a valuable tra cer of extragalactic star formation 
l)Peeters ai!.ll2?)(i^ . However, at this time, not a single 
PA H has been unambiguously identified in the ISM (but 
see lAUamandola efa/l 1)19991) for details of composite 
spectra of the unidentified Infra-red Bands (UIBs) in the 
mid-IR) . At the heart of the problem is the simple fact 
that whil e the AA3.3, 6. 2, 7.7 , 8.6 fc 11.3^m emission 
features l)Peeters etal\ l)2004f )l most likely represent 
excellent tracers of PAH emission, their origins in the 
C-H and C-C stretches, abundant in PAH molecules, 
mean they present little or no diagnostic power. If one is 
to use PAH emission as a key diagnostic on interstellar 
processes, within and outside of our Galaxy, it would 
seem prudent to understand the phenomenon. 

Spatially extended optical emission bands near A5800A 
were first reported by Schmidt et al. ( 1980) but cur- 
rently remain unidentified. A number of detailed descrip- 
tions of the form of the Red Rectangle Bands (RRB) 

^ Hipparchus satellite measurements give a value of 
2.62it2.37mas/yr corresponding to the range ~200-4000pc. 
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ar e available in the l iterat ure. We refer the reader 
to iVan Winckel et all rather than repea t a de- 

tailed description. Following ' Van Winckel et al\ 1)20021) 
we adopt the standard molecular terminology Red De- 
graded (RD) to refer to RRBs which show a sharp 
edge in the blue but an extended tail to the red. The 
A5800A band is a typical RD RRB. Several authors have 
noted the close correspondence of many of the RRBs 
with p rominent DIBs (Scarrott etal. 1992; Sarr e etal\ 
\1995U. rVaii Winckel etalil200'^ . While initial observa- 
tions found strong support for convergence of the RRBs 
to DIBs, the early works were hampered by low resolu- 
tion and low signal-to- noise observations. Later works 
((Glinski & Anderson 2002) indicate that the blue-ward 
shifts of the RRBs, with increasing nebular radius, do 
not appear to asymptotically approach the DIB wave- 
lengths. We concur that the A5800A RRB, in emission 
does not blue shift to the wavelength of the prominent 
A5797A DI B, at least out to a nebu lar radius in excess of 
14.0arcsec. iGlinski fc NuthI (1997^ initially proposed C3 
as the carrier of the A5800A Red Rectangle bands. How- 
ever, more recent observations suggest this is not the 
case ijGlinski fc A ndersonll2002t) . Molecular band struc- 
tures ar c reported towards the centre of the nebula with 
iMivata etali 12004) presenting molecular PAH emission 
in the inner regions of the Red Rectangle in the band 
at ^10fj,m. 

The observations of IGlinski &: AndersonI 1)2002)) . from 
the Densepak IFU on the WIYN telescope, show ra- 
dial evolution of the molecular emission features. How- 
ever, the low surface brightness^ of the Red Rcct- 
angle in the outer region s of the nebula probed by 
IGlinski fc AndersonI l)2002D (10-14arcsec) has hampered 
previous attempts to identify the molecular spectra. 

Recent years have seen considerable progress in lab- 
oratory determination of the spectra of compl ex as- 
trophysically relevant mol ecules (Thaddeus & McQgrJhj 
2001 [lLinnartz eita/.l2004tlSchmidt a/.i2003t.Ding etal. 
2003)) . There has also been progress in the astronomical 
detection of molecular species in the absorption spec- 
traof diffuse clouds l)Maier e^a/.ll2001t lAdamkovics etall 
12003ft and in millimetre-wave emission. However, the low 
surface brightness and complex spatial structure of the 
Red Rectangle have prevented the identification of the 
carrier of its strong molecular emission. 

3. OBSERVATIONS AND DATA REDUCTION 

We used the flames instrument at one of the Nasmyth 
foci of Unit Telescope 2 (UT2 - Kueyen) of the European 
Southern Observatory's Very Large Telescope (ESO vlt) 
located upon Cerro Paranal in northern Chile. The ob- 
servations were undertaken in visitor mode during 29"^ 
and 30*'' December 2004 as part of the Australian guar- 
anteed time awarded for the construction of the OzPoz 
Fibre positioner by the Anglo-Australian Observatory 
dGillingha m . et aL .2003,) . The flames instrument com- 
prises a robotic fibre positioner which provides a range 
of fibre optic feeds to the bench mounted giraffe spec- 
trograph. We used the ARGUS Integral Field Unit (IFU) 
to record spectra of a 2D region of the nebulosity. The 
ARGUS IFU, a rectangular array of 22x14 micro lenses, 
can provide two spatial sampling scales, 0.52arcsec and 

3 HST PC F622W m<17.6mag(Vega) arcsec^ 



0.3arcsec, yielding fields of view of 11.44x7.28arcsec and 
6.6x4.2arcsec respectively. We used the 0.52arcsec lens 
scale for our observations due to the low surface bright- 
ness of the nebula in the regions of primary interest. 

flames is operated in a series of standard user modes. 
We use the HRll mode for our observations, with a 
central wavelength of 5728A and covering a wavelength 
range of 5597-5840A at a dispersion of ~0. 06 A/pixel and 
a resolution element of 2.6pixels (0.156A, R~37,000). 
The detector was a 2kx4k EEV CCD operated in the 
standar d fl ames readout mode. 

Table IXn and Figure fATI detail the observations under- 
taken. Observations were performed during two nights in 
classically scheduled time, 29''^ and SO*'' December 2004. 
While classical mode was not ideal for the observations, 
due to the large range in airmass for the target during 
the night (zenith distance 60°-10°-60°), classical mode is 
a requirement for guaranteed time observations. 

ARGUS attached fiat frames are observed, after stan- 
dard star spectra, twice per night, in order to provide 
relative fibre transmission information and fibre position 
traces. 

Wavelength calibration is derived from standard vlt 
FLAMES daytime calibrations. The giraffe spectro- 
graph offers the ability to observe a number of simul- 
taneous calibration arcs from a Th-Ar lamp during the 
science observations. Unfortunately, a strong line is lo- 
cated close to the primary wavelength region of interest 
leading to significant contamination of the low light level 
nebula spectrum during a 600sec exposures and so simul- 
taneous calibration was not used. 

After discussion with ESO support staff we used the 
ARGUS fast acquisition method whereby acquisition re- 
lies on accurate coordinates for a guide star and will be 
accurate to the order of l-2arcsec. A more accurate, and 
repeatable acquisition would incur a considerable over- 
head for repositioning of guide fibre bundles during the 
night (with each change in PA). The loss of positional 
accuracy has not degraded the observations presented. 

Five base positions were observed within the neb- 
ula. One centred on HD 44179, and four aligned along 
the arms of the nebula, where limb brightening of the 
bi-conical outflow gives the nebula its highest surface 
brightness for a given radial distance, and extending out 
to ~14arcsec. Figure lAll illustrates the IFU locations 
within the nebula. A three point telescope dither was 
performed at each location to allow removal of IFU lens 
artifacts and detector defects, of which there are few in 
the ARGUS/giraffe system. 

ARGUS sky fibres were placed, as required by the AR- 
GUS fast acquisition strategy, in a ring centred on the 
IFU. The ring radius was set to 3arcmins after visual 
inspection of the DSS image of the region showed that 
this would avoid the majority of brighter stars in the 
field regardless of position angle. This avoids the need 
to reposition the fibres between observations which would 
introduce a significant overhead. 

Simultaneous observations with the UVES fibre system 
were not possible due to limits on the proximity of UVES 
fibre placements to the ARGUS IFU. 

3.1. Data reduction 

During the observing period, pipeline data processing 
was not available for the HRll ARGUS setting used. We 
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therefore processed the data using a suite of custom writ- 
ten IDL routines and elements of IRAF. In most cases, 
custom IDL implementations of common iRAF tasks are 
used to allow accurate propagation of error information. 

A 2D bias subtraction was not deemed necessary on 
inspection of the data and hence the data were overscan 
corrected only. On preliminary data reduction it was de- 
cided a known dark current feature seen on the upper 
left corner of the giraffe CCD represented a signifi- 
cant contamination to the low light levels recorded in 
the outer reaches of the nebula and hence subtraction 
of a dark frame was required. A master dark was con- 
structed from 6xl200sec dark exposures. This master 
dark was scaled for the 600sec and ISOOsec exposures as 
no alternate dark frames were available. 

Flat fields are combined using variance weighting with 
outlier rejection. Where three or more science observa- 
tions are available, at a constant dither position, these 
frames are also combined in a similar manner. 

A fibre trace is generated using the iraf apall task, 
applied to the flat field frame. The database trace output 
from the task (a series of S*** order Legendre polynomials) 
was then used to provide the fibre centres for an optimal 
extraction algorithm which accounts for fibre cross talk 
in the spectra. Crosstalk was perceived to be a problem 
during the initial data reduction due to the high surface 
brightness contrast between the inner and outer nebula. 

The extracted fibre flat frames are applied to the sci- 
ence frames after the spectra are extracted in order to re- 
tain the projected relative intensity information between 
the flbres for the extraction process. 

A wavelength solution was derived, using the iraf 
identify and reidentify tasks (fitting a 6*'' order Legendre 
polynomial, with a RMS residual of the order of O.OlA, 
for each fibre) using ThAr arc exposures. The spectra 
are then independently transformed to a common wave- 
length solution. A single 2D transformation, as might be 
used for long slit data, is not possible due to disconti- 
nuities in the wavelength solution as a function of fibre 
number at the boundaries between the sub slits which 
make up the ARGUS slit in the giraffe spectrograph. 

On rectification an estimate of the flat field lamp illu- 
mination spectrum is created from the rectified flat field 
and multiplied back into the spectra to remove the lamp 
signature. 

An absolute flux calibration for the data is not re- 
quired. However, a good relative calibration internal to 
the spectrum is important for molecular modeling. A 
response function is derived via repeated observations of 
the standard star EG-21, a DA white dwarf with few 
spectral features in the wavelength range of interest and 
favourably located on the night sky during observations. 

Spectra are subsequently corrected to a heliocentric 
wavelength solution. 

3.2. Mosaicing the data 

Data from individual cosmic-ray splits (typically 
3xl200sec) are combined using variance weighting and 
outlier rejection. Three dither positions were observed at 
each location in the nebula. In order to align the dither 
data, and to facilitate the ultimate combination of all 
five telescope pointings into a uniform mosaic, we re- 
sample the data onto a common regular Cartesian grid 
after applying the appropriate position angle (PA) ro- 



tation to the IFU internal x/y coordinate system, re- 
trieved from the flames header, to account for the five 
different position angles used during the observations. 
A number of techniques were subsequently investigated 
for aligning individual data cubes. Cross-correlation of 
pseudo-continuum images created from the data cubes 
proved difficult, due to small image overlaps. Ultimately 
alignment was performed using an iterative proce- 
dure, and the iDL imple mentation of the Amoeba algo- 
rithm ijPress a^J |199t|). as applied to a subsection of 
each dispersed data cube. 

The final analysis of the full data set is ongoing. Data 
from only a single arm of the nebula (pointing 2 in Fig- 
ure IA1|) is used to illustrate our model for the A5800A 
emission. Preliminary investigations show that the emis- 
sion in the remaining three nebula arms are qualitatively 
similar. 

4. INTERPRETING THE A5800A BAND. 

A number of origins have been proposed for the 
A5800A emission. Most interpretations propose the 
band profile repres e nts a molecular rotation contour 
Igcarrott etal. 1992: Rouan etal. 1997: Clin ski fc NuthI 
11997) . Glinski fc Anderson (2002) demonstrate that the 
A5800A cannot be generated via Q, R or P branch ro- 
tational emission structures derived from the unknown 
molecular carrier DIB at A5797A. Roiian etal. (19971) 
propose that the RRB has its origins in a supra-thcrmally 
rotating PAH molecule, and go on to argue that a PAH 
of around 40 carbon atoms is implied by the required 
rotational constants. This model assumes that the DIB 
carrier at A5797A is identical to the A5800A RRB car- 
rier and that angular momentum is accumulated due to 
a "rocket effect" induced by photo-dissociation events. 
This effect is required to attain an estimated rotational 
temperature of 450 K for the A5800A carriers. Rotational 
band profiles are always expected of molecular spectra 
since there exist a large number of energetically acces- 
sible rotational states, even when the molecule in ques- 
tion is electronically and vibrationally cold. Thus a sin- 
gle vibronic transition will be derived from a number of 
individual transitions originating in different rotational 
states. Energetically accessible vibrational states, poten- 
tially pumped by photoexcitation, can produce hot-band 
features in an electronic spectrum, where the vibronic 
transitions originate from one of these excited vibrational 
states. Where the potential energy surfaces (PES), gov- 
erning the forces on the nuclei in the electronic states, 
are parallel, the structure of the excitation or emission 
spectrum will be dominated by sequence bands. This 
will be the case for a molecule with an extended chro- 
mophore since the shape of the underlying molecule will 
be largely unaltered upon electronic excitation. 

4.1. Sequence structure 

Based on our study of the radial evolution of the e mis- 
sion structure along the arms of the bi-cone fFigures lAll 
fc IA2p we propose that the A5800A feature arises from 
sequence structure, either associated with the nearby 
A5797A DIB carrier, or an alternate molecule. The ab- 
sence of any perfect correlations between DIBs leads 
one to conclude that if the DIB carriers are molecu- 
lar l)Thorburn e^a^jl2003j) then the Franck-Condon fac- 
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tors for electronic excitation are dominated by the ori- 
gin band. For a molecule in its vibrational and elec- 
tronic ground state, electronic excitation will similarly 
leave the molecule with zero vibration in the electronic 
excited state. The dominance of the origin band in elec- 
tronic excitation is indicative of a transition between 
states with very similar potential energy surfaces. If the 
equilibrium geometry or the vibrational frequencies were 
to change appreciably upon excitation, vibrational pro- 
gressions would be observed. 

If the carrier of the A5800A feature behaves similarly 
to the DIBs, then we would expect any emission from a 
vibrationally excited population of molecules to exhibit 
sequence structure. Sequence structure arises from elec- 
tronic transitions where the molecule remains in vibra- 
tional states of the same character. If the potential en- 
ergy surfaces of two electronic states are parallel, then 
the sequence emission will occur at the same wavelength 
as the origin band. 

In many large molecules, electronic transitions between 
states of increasing vibrational excitation result in tran- 
sitions of progressively lower energy due to the tendency 
for vibrational frequencies to be slightly lower in the ex- 
cited electronic state. In large molecules where electronic 
transitions from the ground state are dominated by the 
origin transition, population of higher vibrational excita- 
tion states results in a red-ward broadening of the emis- 
sion structure. 

This sequence structure interpretation requires the 
molecule responsible for the emission to have an extended 
chromophorc, in order that the electronic transitions do 
not change the structure of the molecule appreciably {i.e. 
the PESs are almost parallel). 

In deference to the current popularity in astronomy of 
PAHs we, in the first instance, follow the trend and pro- 
pose an intermediate size PAH mol ecule as the e mitter . 
Alternative mol ecular fo r ms (C fin - iKroto etoI\ l) 19851) , 
carbon chains - iDouela^ 1)19771) 1. are not considered at 
this time. 

Using the Ultra-High Resolution Facility (UHRF, 
R~600,000 ) of the Anglo- Austrahan telescope (AAT), 
iSarre et all (1995a) show that the A5797A DIB carrier 
fulfils the requirements to generate sequence structure. 
The slightly red-degraded shape of the UHRF spectrum 
shows that the molecule is slightly larger in the excited 
state, indicating a weakening of bonds and thus a prob- 
able slight decrease in the vibrational frequencies of the 
molecule. It is thus suggested that if a vibrationally ex- 
cited population of the carrier of the A5797A DIB were 
to fluoresce (or phosphoresce) it would exhibit sequence 
structure as manifested by red-ward broadening of the 
DIB band shape. However, as the A5800A feature does 
not asymptotically approach the A5797A DIB position, 
the simplest assumption to be made is that they arise 
from a different carrier. For the bands to arise from the 
same carrier we require that the coldest emission be com- 
pletely self-absorbed by the nebula. 

4.2. PAH Model 

In order to explore the sequence band hypothesis, we 
performed normal mode analyses on PAH cations of var- 
ious sizes in ground and electronically excited states. 
Cations, rather than neutrals, were chosen for these 
quantum calculations so that they exhibited strong tran- 



sitions in the visible region ijBrechignac fc PinollT999(l . 
Ground state geometries and normal mode frequencies 
of the naphthalene, pyrene, tetracene, perylene and 
coronene cations were obtained with the B3LYP density 
functional and 6-31G basis set using the GAUSSIAN 03 
(,Frisch et al. 2003) suite of programs. Time-dependent 
Density Functional Theory (DFT) calculations were per- 
formed on the optimised ground state structures to yield 
excited states. The electronic transition in the visible 
region with the largest oscillator strength was chosen. In 
order to replicate this excited state, the ground state or- 
bital occupancy was altered to resemble the leading con- 
figuration of the selected excited state. Geometry opti- 
misations and normal mode analyses were performed on 
these excited states for naphthalene, pyrene, tetracene 
and perylene. 

The vibrational frequencies in the ground {v'/} and 
excited states {i/-} were used to generate transition en- 
ergies with the origin transition (arbitrarily) centred at 
A5800A, in order that we be able to compare the band 
shape with the new RRB observations. At the time of 
writing, only indicative origin positions may be calcu- 
lated by quantum theory for molecules of the size being 
invoked here. By way of bench marking, pentacene is seen 
to absorb at A5362A llHeineckell998f) but is calculated by 
TD-DFT to absorb at 6190A (Reilly & Schmidt, private 
communication). Sequence band profiles were calculated 
in the manner laid out in appendix IXI 

For an example of a n observed mo l ecular sequence 
structure see Figure 2 of lLinnartz etal\ l)2004D . 

4.3. DFT Modeled emission spectra 

It must be emphasised that the modelling of emis- 
sion spectra is, in essence, an exploration of the pa- 
rameter space available to the molecules responsible for 
the A5800A emission. As such, by invoking naphthalene, 
pyrene, tetracene and perylene cations we are exploring 
the expected vibrational sequence structure of a molecule 
of that size and shape, with an extended chromophore, 
rather than the species itself. 

Notwithstanding the fact that naphthalene cation is 
likely photodestroyed in harsh circumstellar environ- 
ments, it serves as a starting point for the exploration 
of our parameter space. As can be seen in figure IX3l the 
modeled naphthalene spectra are far too structured and 
cover too large a wavelength region to be responsible for 
the RRBs. What's more, the density of states does not 
increase fast enough with energy to bring about a shift 
in the maximum of the emission spectrum. 

On moving to pyrene, we see that the spectra become 
less resolved due to an increased number of emissive 
states. At first glance it would appear that pyrene be- 
haves in a very similar manner to the new RR spectra. 
However, the shift in the emission peak for pyrene is too 
large to be considered a se rious candidate for the RR 
emission carrier (c/. Figures IX2l fc IX3|) . 

Tetracene is calculated to have small vibrational fre- 
quency shifts as compared to naphthalene and pyrene. 
As such it yields narrow sequence structure features 
which are broadly in keeping with the behaviour observed 
in the RR. 

Perylene behaves in a similar way to pyrene yet its 
sequence structure is calculated to be somewhat domi- 
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nated by a small frequency mode with a relatively large 
shift upon excitation. The greater number of Boltzmann 
accessible vibrational states accessed in the electronic ex- 
cited state brings about a shift in the maximum of the 
emission spectrum to longer wavelengths. However, even 
for perylene the magnitude of the red-shift of the calcu- 
lated sequence structure is far greater than that observed 
in the RRBs. In order that the RRBs be explained by 
sequence structure, we require a larger molecule with a 
much larger density of states. 

The size of the molecule required is such that the mod- 
eling of its sequence structure explicitly using DFT be- 
comes intractable. While it may be possible to obtain 
ground state electronic structures and vibrational fre- 
quencies of very large molecules, the excited state elec- 
tronic structures become increasingly complicated, being 
brought about by mixing several excitations out of the 
ground state. As such, representation of an excited state 
with an altered orbital occupancy is unreliable. Optimiz- 
ing excited states by using TD-DFT is beyond the scope 
of contemporary computer codes for the size of molecules 
being employed here {vide infra). 

For molecules the size of coronene and larger, excited 
state frequencies were chosen guided by the behaviour of 
the smaller systems at the DFT level of theory. Figure 
IA4l demonstrates that most frequency shifts upon excita- 
tion fall within an envelope (shaded). Randomly select- 
ing shifts from this envelope changes the nature of the 
sequence structure very little, as shown in figure IX5l for 
pyrene. Similarly, calculating the ground state frequen- 
cies using molecular mechanics (MM) produces rather 
similar results. The correlation between MM a nd D FT 
calculated frequencies is tight, as shown in Figure lXHl To 
investigate systems larger than coronene we have invoked 
a purely molecular mechanical model, and assumed that 
these structures behave in the same manner observed by 
DFT for the smaller PAH cations. Coronene frequen- 
cies were calculated with both DFT and MM, with both 
results presented. 

5. MOLECULAR MECHANICAL MODEL 

Using parameter s taken from the CHARMM force-field 
l)Brooks etaUW Sl) , vibrational frequencies for electronic 
ground states of large PAHs were generated by diagonal- 
ization of the mass-weighted Hessian matrix. Excited 
state vibrational frequencies were chosen from the enve- 
lope of DFT shifts as outhned above. Sequence structure 
was then simulated in an identical manner to the implicit 
DFT results. 

5.1. Molecular Mechanical model results 

Using a molecular mechanical (MM) model, we ex- 
plore a range of larger PAH molecules, including pyrene, 
coronene, hexabenzocoronene, circumcoronene and dicir- 
cumcoronene. 

Figure IA5I compares the modeled emission spectra of 
pyrene at 200 K from the DFT model and the MM model. 
The effect of increasing the size of the molecule in the 
MM model is that the density of states increases very 
rapidly with energy and thus the sequence structure is 
shifted to the red. Spectra for the coronene, hexaben- 
zocoronene, circumcoronene and dicircumcoronene at a 
range of temperatures are shown in Figure 1X71 



Given the similarity of the sequence structure pre- 
dicted by the MM model it is thus difficult to pin-down 
the exact size of the A5800A carrier. 

Despite the molecules invoked here displaying similar 
red-shifted emission spectra, the red-shift of the emission 
maximum was found to be very temperature sensitive. 
As such the MM modeled temperature of the Red Rect- 
angle due to sequence emission is relatively insensitive 
to carrier(!). We found that the broadest emission spec- 
tra observed in the Red Rectangle could be explained in 
terms of molecules at vibrational temperatures of ~90K. 
Modeled MM spectra for our fiducial molecule, dicircum- 
coronene, are superimposed on A5800A RRB in Figure 
IA8I An initial estimate of the radial vibrational tem- 
perature of the A5800A band carrier is given in Figure 

5.2. Deficiencies of the MM model 

While the redward shifts in the maxima of emission 
are reproduced with little trouble, our model may neces- 
sarily over-estimate the size of the PAH required for the 
observed (putative) sequence structure. This is because 
we select the range of frequency shifts from DFT calcula- 
tions necessarily performed on small molecules. One may 
expect that larger molecules will exhibit smaller shifts 
upon excitation as a consequence of the electronic tran- 
sition being more delocalized. The effect of this will be to 
pull the sequence band contour into the origin band po- 
sition, reducing the amount of red shifting. As such, the 
observed structure in the RR may be due to a smaller 
molecule than dicircumcoronene, at a higher tempera- 
ture. 

Additionally, our model only contains a single temper- 
ature population at this time. The effect of a compound 
population along the line of sight, with a restricted range 
temperature, would be to smear out the emission peaks 
somewhat. 

5.3. Should sequence structure be observed? 

iSarre et'al\ l)1995bD note that one might expect se- 
quence structure to be a feature of molecular emission 
bands. If one accepts that large PAH molecules exist 
in appreciable quantities in the Red Rectangle, and that 
they will be undergoing collisions and photoexcitation, 
then one concludes that they will exhibit a significant 
vibrational temperature. While many vibrational modes 
can radiate energy (indeed, the 3.3/im and 7.6/im UIB 
emission is proposed to originate from PAH C-H and C-C 
stretches) there will be many modes which cannot cou- 
ple to electromagnetic radiation. As such, vibrational 
energy can be "stored" in the PAH for an appreciable 
time. This phenomenon and the high photon flux in 
the nebula will bring about a vibrational temperature. 
This temperature may exceed 100 K, and electronically 
excited molecules will radiate to vibrational wavefunc- 
tions with nearly identical character to that in the ex- 
cited state. We conclude that a large PAH molecule in 
an environment such as the RR will inevitably exhibit 
sequence structure in its emission. 

5.4. De-hydrogenation 

Kokkin & Schmidt (in prep.) demonstrate that the ef- 
fect of de-hydrogenation on the in-plane electronic tran- 
sitions in large PAH molecules is primarily to shift the 
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band a small amount (of the order of a few ~10A) al- 
though bond-length alteration induced by dehydrogena- 
tion may be important in systems as large as decacyl. 
Since in-plane transitions are the strongest and only in- 
volve molecular orbitals of 7r-symmetry {A"), dehydro- 
genation has little effect on the position of the excitation. 
Out-of-plane transitions of PAH radicals are expected to 
be significantly shifted from the strong in-plane transi- 
tions of the neutral close-shell parent, and such spectra 
will also exhibit a degree of vibrational structure due to 
the localization of the lone-pair chromophore. 

If the A5800A band should indeed be attributed to 
emission from a 45-100 carbon atom PAH (such as di- 
circumcoronene) then one may postulate a second band, 
at longer wavelengths, derived from the de-hydrogenated 
molecule. Such bands would almost certainly be ob- 
served in vibrationally hotter states since one expects 
these bands to be more prominent closer to the central 
star, an energy source capable of removing the hydro- 
gen atom from the PAH. They would not be observed as 
sharp bands such as the RD e mission feature. Suc h bands 
would better resemble the I Van Winckel etal\ l)2002f) 
Symmetric (S) bands. The Red Rectangle spectrum is 
replete with such features. Th us we tentatively propos e 
the second feature, at A5826.5 (' Van Winckel etall\200^ . 
the extreme red limit of our ARGUS spectra, as sequence 
emission from a dehydrogenated PAH responsible for the 
A5800A band. 

5.5. Other Red Degraded (RD) bands. 

If one accepts a sequence structure explanation for the 
A5800A RRB, it is tempting to propose similar expla- 
nations for the remaining RD bands in the Red Rectan- 
gle. A number of RD bands, each of which can plausibly 
be assigned an associated S band de rived from the de- 
hydro genated molecule, are listed bv lVan Winckel etall 
l)2002j) . Such bands would not arise in a common carrier 
but would instead arise in similar 45-100 carbon atom 
PAH molecules. Where a PAH has more than one sym- 
metry unique dehydrogenation site, a single RD band 
could have associated with it more than one S band. 

5.6. Symmetry and candidate selection 

When suggesting candidate carriers for the RRBs, 
one is naturally tempted to choose highly symmetric 
molecules in order to avoid the plurality associated with 



arbitrary choices. However, there are other good rea- 
sons to suggest a highly symmetric carrier. The Born- 
Oppenheimer breakdown which prevents many large 
molecules from fluorescing is facilitated by a large back- 
ground of vibronic states of the same symmetry. By in- 
voking a more symmetric candidate, the number of back- 
ground states of the same symmetry naturally falls, as 
there are more irreducible representations of the molec- 
ular point group. As such, it is more likely that a highly 
symmetric molecule will fluoresce than a less symmetric 
one. 

6. CONCLUSIONS 

We have demonstrated that the radial variation of the 
profile of the A5800A Red Rectangle band is well mod- 
eled as molecular sequence structure arising in a PAH 
molecule with 45-100 carbon atoms. Our model has the 
attractive properties of : 

• Not requiring an exotic molecule, once one accepts 
the presence of PAHs in the ISM. 

• Not requiring high nebular temperatures. 

• Not requiring an exotic emission mechanism. 
While sequence structure does not appear promi- 
nently in the astronomical literature, it is a simple 
molecular emission phenomenon. 

If sequence structure is the origin of the A5800A RRB 
then it is natural to extend the hypothesis to the other 
Red Degraded RRBs, invoking sequence structure for ad- 
ditional molecules. 
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APPENDIX 

DFT CALCULATIONS OF BAND PROFILES 
Assuming no change in vibrational state, the photon energy for sequence transition j, Tj, is simply 

3N-6 
i=l 

where To is the energy of the origin transition, N is the number of atoms in the PAH cation and is the number 
of quanta of vibrational energy in mode i. Simulated emission spectra were produced by convolving the predicted 
transitions energies with a Gaussian profile simulating the rotational structure of each sequence band. 

The vibrational population in the excited electronic state was assumed to have an exponential dependence on 
vibrational energy of the ground vibrational state j 



P, = exp(-/3£;;), 



(A2) 
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where 

3iV-6 

E;^hJ2 "y^r- (A3) 

i=l 

The simulated emission spectrum /(A) is thus a sum over sets of quanta (vibrational states) 

J(A) = Pj exp(-a(A - hc/T^f) (A4) 

3 

where a controls the FWHM of the gaussian (rotational) profile. 

In cases where the number of states needed to converge the shape of the calculated sequence emission profile was 
prohibitively large, a Monte Carlo sampling technique was employed. Here, states, j, were chosen as a Monte Carlo 
Markov Chain and added to the simulation one at a time (starting at the ground vibrational level). The probability 
of accepting a step, in the space of the number of quanta in each vibrational mode, was given by the Metropolis 
technique, whereby all steps down in vibrational energy were chosen with unit probability and those up in energy were 
accepted with probability Pj — exp(— /JE""). Since the probability of state j being occupied is naturally taken into the 
selection algorithm, the emission spectrum was taken as a simple sum over the sampled states. 

/(A) = cxp(-a(A - hc/Tjf). (A5) 

3 

Typically, 10000 states were enough to converge the shape of the emission band profile. 
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TABLE Al 

Observation log, 20041229-30, VLT-UT2 ARGUS HRll. 



Pointing 


RA/Dec (J2000) 


Pos. Ang. 


Exp. time 


Date 


PO 


06:19:58.14 -10:38:14.0 


12 


3x2x600sec 


2005-12-29 


PI 


06:19:58.32 -10:38:23.3 


165 


3x2xl800sec 


2005-12-29 


P2 


06:19:57.68 -10:38:20.7 


230 


2x3xl200sec 


2005-12-29 








Ix3xl200sec 


2005-12-30 


P3 


06:19:57.95 -10:38:04.5 


340 


3x3xl200sec 


2005-12-30 


P4 


06:19:58.63 -10:38:06.4 


43 


3x3xl200sec 


2005-12-30 
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Fig. A1. — The l ocations and orientations of base telescope pointings and the ARGUS IFU are shown overlayed on the HST PC image, 
at F622W, of Cohen etal. 12004 ). The large cir cles are drawn at a radius of 10.7 and 14.3arcsec, which correspond to the fiducial WIYN 
Densepak spectra of lGlinski fc AndersonI 120021) . 
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5770 5780 5790 5800 5810 5820 5830 5770 5780 5790 5800 5810 5820 5830 

Wavelength(A) Wavelength(A) 
Fig. A2. — As an example of the strong evolutio n of t he A5800A band, as a function of radial distance from the central star, a sequence 
of spectra from the south western arm (P2 in Figure lXTl of the nebula is shown. The radial distance from HD 44179 is indicated to the left. 
The spectra are vertically baseline offset for clarity. Spectra in the left-hand Figure have been continuum subtracted, while the identical 
spectra in the right-hand Figure have been normalised by the continuum. In this illustration the continuum has been estimated as a simple 
constant value over the spectral range AA5780-5790A. A vertical dashed line illustrates the location of the A5797A DIB. Tick marks indicate 
the peak of the A5800A band within each spectrum. 



1 200 K [COt I 300 K 



300 K 




[c ccor 1^ [Q ccor 200 k [oxo]^ 

ri . looK n ., 200K n , 300 k 

5720 5760 5800 5840 5880 5720 5760 5800 5840 5880 5720 5760 5800 5840 5880 

Wavelength (A) Wavelength (A) Wavelength (A) 

Fig. A3. — Calculating the sequence structure emission profile using Density Functional Theory (DFT), we conclude that small PAHs 
alter their structure considerably on excitation. Consequently, the red-shift in the band peak is too great and too much structure is resolved 
within the band profile for molecules of this size to represent the A5800A band evolution. 
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Fig. A4. — A plot of ground state DFT vibrational frequencies (t/") against their respective shifts in the excited electronic state (Ai/). 
For molecules larger than perylene explicit DFT calculations for the ground and excited states are intractable. In order to model sequence 
structure emission for larger molecules, to demonstrate how such emission would be observed, we assume a series of randomly realised 
frequency shifts drawn from the envelope of parameter space populated by a number of PAH molecules for which explicit DFT calculations 
of both the ground and excited states were possible. 
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Fig. A5. — For the PAH pyrene explicit DFT calculations of the ground and excited state spectra are possible. Here we contrast the 
explicit DFT calculations (DFTe) with the implicit model outlined in the text using the ground states frequencies derived from the DFT 
calculations (DFT^) and the MM model (MM^). 
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Fig. A6. — For molecules larger than coronene DFT calculations failed to converge on acceptable ground state frequencies and we are 
forced to resort to a Molecular Mechanical model to generate the ground state frequencies. We confirm the validity of the MM model 
by comparing the resulting ground states to those found via DFT for small PAIf molecules for where DFT is possible. A one-to-one 
correspondence is indicated by the dashed line. The large frequencey gap between 1700 cm"'^ and 3100 cm~^ is due to the mass difference 
between carbon and hydrogen, the region around 3100 cm~^ being due to hydrogen stretching motions. 
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Fig. A7. — Sequence structure, as a function of temperature, is calculated for a number of large PAH molecules using the MM model. 
For coronene, the largest molecule for which we are able to determine reliable frequencies for the ground state via DFT , the implicit DFT 
model (DFTi) is shown by the shaded region. Sequence structure, exhibited as a natural consequence of the molecular structure of large 
PAH molecules, closely resembles the radial variation in the spectra of the red degraded (RD) Red Rectangle bands such as the A5800A 
band. 
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Fig. A8. — A comparison of the MM spectra, for the molecule dicircumcoronene, with the astronomical spectra is shown. Choosing 
dicircumcoronene as our fiducial molecule, the only variable in the model is the temperature governing the distribution of excited states. 
Sequence structure closely mimics the behaviour of the observed spectra. We believe the feature to the right of the main band to be a 
related but large ly independent feature of the spectrum which is not to be explained by sequence structure from the ground state transition 
(see section [5.41 . 
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Fig. A9. — Using the Molecular Mechanical model developed in the text and our fiducial molecule, dicircumcoronene, one can us e the 
peak shift for the sequence bands to estimate a radial temperature variation for the nebula, in this instance along arm 2 from Figure IXTl 



